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ABSTRACT 

Context. To add to the growing collection of well-studied double periodic variables (DPVs) we have carried out the first spectroscopic 
and photometric analysis of the eclipsing binary DQ Velorum to obtain its main physical stellar and orbital parameters. 
Aims. Combining spectroscopic and photometric observations that cover several orbital cycles allows us to estimate the stellar prop- 
erties of the binary components and the orbital parameters. We also searched for circumstellar material around the more massive star. 
Methods. We separated DQ Velorum composite spectra and measured radial velocities with an iterative method for double spectro- 
scopic binaries. We obtained the radial velocity curves and calculated the spectroscopic mass ratio. We compared our single-lined 
spectra with a grid of synthetic spectra and estimated the temperature of the stars. We modeled the V-band light curve with a fitting 
method based on the simplex algorithm, which includes an accretion disc. To constrain the main stellar parameters we fixed the mass 
ratio and donor temperature to the values obtained by our spectroscopic analysis. 

Results. We obtain a spectroscopic mass ratio q = 0.31 ± 0.03 together with donor and gainer masses Ma = 2.2 ± 0.2 M a , 
M g = 7.3 ± 0.3 M , the radii R d = 8.4 ± 0.2 R e , R g = 3.6 ± 0.2 R and temperatures T A = 9400 ± 100 K, T g = 18 500 ± 500 K for 
the stellar components. We find that DQ Vel is a semi-detached system consisting of a B3V gainer and an A1III donor star plus an 
extended accretion disc around the gainer. The disc is filling 89% of the gainer Roche lobe with a temperature of 6580 ± 300 K at 
the outer radius. It has a concave shape that is thicker at its edge (d c = 0.6 ±0.1 R ) than at its centre (d c = 0.3 ±0.1 R Q ). We find a 
significant sub-orbital frequency of 0. 19 d~ l in the residuals of the V-band light curve, which we interpret as a pulsation of an slowly 
pulsating B-type (SPB) of a gainer star. We also estimate the distance to the binary {d ~ 3.1 kpc) using the absolute radii, apparent 
magnitudes, and effective temperatures of the components found in our study. 

Key words. Stars: binaries, stars: eclipsing, stars: fundamental parameters, stars: early-type 

1. Introduction also used low-resolution spectra obtained at the Cerro Tololo 

„. „ _ „ . „ T , , ,__ „. _ inter-american observatory with the image-tube spectrograph at 

£ ' ^TTt^^TlZ fi^' DQ/elorum (TYC 8175- ^ Yale 1 . m ^ t ? an estim g ate of J spec |/ type 

3 ■ 333-1, ASAS J093034-501 .9) is a Galactic echpsmg binary afj A3y Howeyer> n P o wafj made tQ ^ 

or Algol-type with an orbital period or 6.08337 d (Horfmeister velocities 

E94|. Considering that DQ Vel is a bright system (V~ 10.7), it More recentl Michalska et al .rjoij discovered an additional 

is still poorly analysed in detail long-term variability of 188.9 d using ASAS photometric 

st , u , d y of *l|y sten ^ U ' lng Photographic plates was made by data Thig lacefj ^ em ^ ^ flew of yaj iaWe star§ 

van Houten ([1953, who determined the following ephemens for called doMe pmW£c variables (hereafter DPVs) . This group 

the primary minimum. of interacting binaries is characterised by two photometric 

T min = JD 2430881.20(3) + 6.08337(13) x E . (1) variabilities linked one to each other with a period ratio of 

about 33 (Mennickent et al. 2003 , Mennickent & Kolaczkowski 

Milone (1986), using UVB photometry, indicated that DQ Vel 1201 Oi l. The short-term variability corresponds to the orbital 

shows a /? Lyrae-like light curve and in vestigated the O'Connell motion of the binary while the long-term variability is still 

effect on the system (0'Connell ll951l l. studying the differential not completely understood. DPVs have been interpreted as 

light between maxima of the light curves. He concluded that semi-detached binaries showing cycles of mass loss into the 
the O'Connell effect is positive and stronger in the U band. He 
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interstellar medium (Mennickent et al. 20081 120 12bt . 
To study these interesting objects, an intense observational 
effort is being carried out to find DPVs both in the Galaxy and 
the Magellanic Clouds using the microlensing surveys OGLE 
and MACHO (Mennickent et al. 120051 Poleski et al. [20T0b , as 
well as in the ASAS project. At present, there are more than 
150 DPVs discovered in the Magellanic Clouds and 1 1 systems 
(including DQ Vel) in the Milky Way (Mennickent et al. l2012al) . 
Results of multiwavelength spectroscopy and/or photometry 
have made it possible to study DPVs in more detail, such as AU 
Mon (Desmet et al. 120101 Dju rasevic et al. l20T0b , V393 Sco 
(Menn ickent et al. 120 101 1201 lal . and LP Ara (Mennickent et al. 
1201 lbl) . 

Here, we present the first optical spectroscopic analysis of DQ 
Vel using 46 high-resolution spectra that cover several orbital 
cycles. Our analysis includes the separated components with a 
determination of their radial velocity curves, an estimation of 
the spectroscopic mass ratio, stellar and system parameters, and 
the spectral type of the two components. 

To constrain the main stellar parameters we also include in our 
analysis V-band light curves obtained from the public ASAS 
archive and VIJK photometric data obtained by Michalska et al. 
(120101 ) at the REM 0.6-m telescope at La Silla. 



2. Observations 

To investigate the fundamental properties of DQ Vel we col- 
lected a series of high-resolution optical spectra between April 
2008-2011. Most of the spectra were taken with the 1.2m 
Leonard Euler Swiss telescope at the La Silla Observatory using 
the CORALIE echelle fiber-fed spectrograph. In our analysis we 
also included five spectra obtained in June 2008 and May 2009 
from the 2.5m Irenee du Pont telescope echelle spectrograph at 
the Las Campanas Observatory as well as two spectra taken with 
the echelle spectrograph FEROS mounted at the ESO/MPI 2.2m 
telescope at La Silla. A detailed summary of our spectroscopic 
observations is given in Table 1 . The orbital phases <p„ were cal- 
culated using the ephemeris obtained by Michalska et al. (120101) : 
they are detailed in Section 3.3. The echelle spectrograph at Du 
Pont has a wavelength coverage from ~ 3700 to 7000 A. The in- 
strument resolution is 45 000 for one arc-second slit. These data 
were reduced with the standard procedures of de-biasing, flat- 
fielding, background subtraction, extraction of one-dimensional 
spectra, wavelength calibration, cosmic spike removal and con- 
tinuum normalisation using the Image Reduction and Analysis 
Facility IRAF Q. The reduced spectra were then subjected to 
barycentric corrections. 

The high-resolution FEROS spectrograph (R~48 000) has a 
slightly broader spectral range than CORALIE, from 3500 to 
9200 A (Kaufer et al. [19991 120001 The FEROS data presented 
here were previously reduced by the FEROS Data Reduction 
System (DRS), which includes barycentric corrections and the 
standard procedures mentioned above except for continuum nor- 
malisation, which was made using the continuum task of IRAF. 
The CORALIE spectrograph was developed to measure pre- 
cise stellar radial velocities. The maximum spectral resolution 
of about 100000. The instrument has an automatic reduction 
software (Baranne et al. 1 19961 1 that provides fully calibrated 

1 IRAF (http://iraf.noao.edu/) is distributed by the National Optical 
Astronomy Observatories and operated by the Association of 
Universities for Research in Astronomy Inc., under cooperative agree- 
ment with the National Science Foundation. 



Table 2. Detailed photometric data used in this paper. First col- 
umn is the observing window and /V is the number of frames. 



BJD-2450000 Observatory Telescope Filter N 



1998.8357-2885.4690 


La Silla 


REM 


V 


519 


1998.8372-2899.4789 


La Silla 


REM 


I 


532 


1998.8360-2865.5494 


La Silla 


REM 


J 


211 


2000.8472-2855.4982 


La Silla 


REM 


K 


135 


1868.8363-4971.6065 


Las Campanas 


ASAS 


V 


583 



images and computes barycentric velocities. Just as with the 
FEROS spectrograph, CORALIE uses a simultaneous Thorium 
lamp spectrum during the science exposures, which allows a reg- 
ular wavelength calibration. Here we used the CORALIE spec- 
tra reduced with this automatic pipeline. Cosmic spikes were re- 
moved manually. It should be noted here that we found a small- 
amplitude (lower than 1% of the total flux) continuum modu- 
lation, most likely produced by an imperfect automatic extrac- 
tion of the orders. To attenuate or remove this modulation we 
subdivided each CORALIE spectrum into different wavelength 
intervals. For each interval we applied the continuum task us- 
ing the spline3 function with a corresponding order (different 
in each interval), to create a "continuum normalised "interval. 
We repeated this procedure for all spectra using the same wave- 
length intervals, same function, and same orders as used for 
the first spectrum. Finally, we combined the corresponding in- 
tervals of each spectrum to create a final normalised spectrum. 
Although the data quality significantly improved, there still re- 
main some residuals, which were further attenuated by applying 
a high-frequency Fourier filter to the spectra during the cross- 
correlation process, as will be explained in Section 3.1. 
Our photometric analysis includes a V-band light curve obtained 
from the public ASAS database 0. The observations include 
583 frames taken in different epochs during November 2000- 
2009. Typical errors on these data are less than 0.06 mag. The 
VIJK photometric data used here were obtained with the REM 
0.6m telescope at La Silla, which were previously reported by 
Michalska et al. (120101) . who decoupled the short- and long-term 
variabilities using the Fourier decomposition technique. The ob- 
servational campaign included 1051 frames in the V-I bands and 
345 frames in the J-K bands. Details of the photometric data are 
presented in Table 2. 



3. Results 

3.1. Spectral separation and radial velocity curves 

To investigate if DQ Vel is a single or double line spectro- 
scopic binary we analysed the 46 spectra listed in Table 1. 
We detected known absorption lines and calculated their cor- 
responding Doppler shifts. Excluding the intense and broader 
Balmer lines that seem to contain several components, we 
selected the stronger and unblended lines: Hel 4387.929 A, 
MgII4481.126A, Fel 4957.596 A, Hel 5875. 6 14 A, the silicon 
doblet Sill 6347 . 1 03 A/637 1 .359 A, and Hel 6678. 1 5 1 A. We di- 
rectly measured each line position fitting a Gaussian profile us- 
ing the IRAF task splot and calculated their corresponding radial 
velocities. We measured each line at least four times to obtain an 
estimate of the error. We found that all selected Hel lines seem 
to belong to the more massive star (hereafter gainer) while the 
metal lines apparently belong to the less massive star (hereafter 

2 www.astrouw.edu.pl/asas/ 
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Table 1. Detailed optical spectroscopic observations of DQ Vel. </> is the orbital phase calculated using the Michalska et al. (1201 Oi l 
ephemeris (see section 3.3). The signal-to-noise ratio was measured using the 5000-5800 A range at the continuum level. 



Observatory 


Telescope/Instrument 


UT date 


BJD-2450000 


<b„ 

TO 


Exp time(s) 


S/N 


La Silla 


Euler/CORALIE 


2008-04-05 


4562.5807 


0.8061 


2500 


53 


La Silla 


Euler/CORALIE 


2008-04-05 


4562.6105 


0.8110 


2500 


57 


La Silla 


Euler/CORALIE 


2008-04-05 


4562.6403 


0.8159 


2500 


55 


La Silla 


Euler/CORALIE 


2008-04-07 


4564.6168 


0.1408 


2500 


60 


La Silla 


Euler/CORALIE 


2008-04-07 


4564.6466 


0.1457 


2500 


55 


La Silla 


Euler/CORALIE 


2008-04-08 


4564.6764 


0.1506 


2500 


56 


La Silla 


Euler/CORALIE 


2008-04-08 


4565.5829 


0.2996 


2500 


63 


La Silla 


Euler/CORALIE 


2008-04-08 


4565.6127 


0.3045 


2500 


63 


La Silla 


Euler/CORALIE 


2008-04-08 


4565.6425 


0.3094 


2500 


67 


La Silla 


Euler/CORALIE 


2008-05-22 


4609 545 1 


0.5263 


2500 


58 


La Silla 


Euler/CORALIE 


2008-05-22 


4609.5750 


0.5312 


2500 


64 


La Silla 


Euler/CORALIE 


2008-05-22 


4609 6258 


0.5395 


2500 


57 


La Silla 


Euler/CORALIE 


2008-05-23 


4610.5343 


0.6889 


2500 


57 


La Silla 


Euler/CORALIE 


2008-05-23 


4610 5641 


0.6938 


2500 


58 


La Silla 


Euler/CORALIE 


2008-05-23 


4610.5940 


0.6987 


2500 


55 


La Silla 


Euler/CORALIE 


2008-05-24 


4611.5221 


0.8512 


2500 


49 


La Silla 


Euler/CORALIE 


2008-05-24 


4611.6078 


0.8653 


2500 


46 


La Silla 


Euler/CORALIE 


2008-05-24 


4611.6377 


0.8702 


2500 


42 


La Silla 


Euler/CORALIE 


2008-05-25 


4611.6675 


0.8751 


2500 


61 


La Silla 


Euler/CORALIE 


2008-05-25 


4612.5771 


0.0247 


2500 


34 


La Silla 


Euler/CORALIE 


2008-05-25 


4612.6069 


0.0296 


2500 


32 


La Silla 


Euler/CORALIE 


2008-05-25 


4612.6367 


0.0345 


2500 


31 


La Silla 


Euler/CORALIE 


2008-12-19 


4819.7745 


0.0843 


1500 


42 


La Silla 


Euler/CORALIE 


2008-12-22 


4822.7697 


0.5767 


1500 


48 


La Silla 


Euler/CORALIE 


2009-04-15 


4936.6295 


0.2932 


2500 


62 


La Silla 


Euler/CORALIE 


2009-04-15 


4936.6594 


0.2981 


2500 


60 


La Silla 


Euler/CORALIE 


2009-04-15 


4936.6892 


0.3030 


2500 


59 


La Silla 


Euler/CORALIE 


2009-04-16 


4937.6087 


0.4542 


2500 


49 


La Silla 


Euler/CORALIE 


2009-04-16 


4937.6386 


0.4591 


2500 


55 


La Silla 


Euler/CORALIE 


2009-04-16 


4937.6684 


0.4640 


2500 


53 


La Silla 


Euler/CORALIE 


2009-04-17 


4938 6259 


0.6214 


2500 


57 


La Silla 


Euler/CORALIE 


2009-04-17 


4938.6701 


0.6287 


2500 


54 


La Silla 


Euler/CORALIE 


2009-04-17 


4938.6999 


0.6336 


2500 


49 


La Silla 


Euler/CORALIE 


2009-05-17 


4968.5220 


0.5358 


2400 


50 


La Silla 


Euler/CORALIE 


2009-05-17 


4968.5507 


0.5405 


2400 


47 


La Silla 


Euler/CORALIE 


2009-05-17 


4969.5085 


0.6980 


2500 


53 


La Silla 


Euler/CORALIE 


2009-05-18 


4969.5378 


0.7028 


2400 


55 


La Silla 


Euler/CORALIE 


2009-05-18 


4970.4967 


0.8604 


2400 


49 


La Silla 


Euler/CORALIE 


2009-05-19 


4970.5253 


0.8651 


2400 


44 


La Silla 


MPI/FEROS 


2010-03-28 


5283.5531 


0.3214 


1500 


67 


La Silla 


MPI/FEROS 


2010-03-30 


5285.5668 


0.6524 


1500 


62 


Las Campanas 


Du Pont/echelle 


2008-06-13 


4631.4735 


0.1309 


900 


27 


Las Campanas 


Du Pont/echelle 


2008-06-13 


4631.4854 


0.1329 


900 


46 


Las Campanas 


Du Pont/echelle 


2009-05-14 


4966.4639 


0.1975 


900 


60 


Las Campanas 


Du Pont/echelle 


2009-05-15 


4966.5077 


0.2047 


900 


49 


Las Campanas 


Du Pont/echelle 


2009-05-16 


4967.5020 


0.3681 


1200 


60 



donor). Averaging the Doppler shifts measured for each line, we 
calculated mean radial velocities (RVs) for the gainer and donor 
star in each spectrum. The resulting RV measurements are plot- 
ted in Figure 1 . Preliminar inspection of these RV curves shows a 
defined sinusoidal curve for the donor star. The orbit seems to be 
circular; if eccentricity is present, its value is negligible within 
the errors. 

We observed more scatter on the gainer velocities. We attributed 
this to the fact that the helium lines are broader than the donor's 
metal lines, which results in a higher uncertainty on the deter- 
mination of the line centre. Possible causes for the broadening 
of the helium lines and/or for the large scatter observed on the 
gainer RVs are either a hot rapidly rotating star, as expected 
for a gainer that is accreting angular momentum together with 
the material, or line profiles blended by the contribution of a 
non-symmetrical pseudo-photosphere or disc. Yet another cause 



could be that the gainer is a pulsating star (See section 3.5). 
Using the last calculated RVs, we decided to separate the in- 
dividual spectra to obtain more precise RVs. We used a simple 
iterative method originally proposed by Marchenko et al. ( 1 998), 
which was used more recently by Gonzalez & Levato (2006), 
who applied this method for the orbital analysis of the binary 
HD 1435 1 1 . In general terms, this method is able to separate both 
components, providing two higher signal-to-noise (S/N) single- 
line spectra templates, which can be used to calculate individual 
RVs, for instance through a cross-correlation process (Griffin 
119671 Hill I1993I ). The Fourier cross-correlation process (here- 
after CCR) convolves two functions (spectra against the tem- 
plate) using a Fourier transform algorithm to find the so-called 
cross-correlation function ( CCF). Considering that ~ 80% of our 
data are from the CORALIE spectrograph and that the accuracy 
of a relative velocity measured with the CCR process increases 
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300 



200 



100 



-100 

-200 
0.0 



X MgIM481/FeI4957/SiII6347/6371 
• HeI4387/5875/6678 
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\ 



0.2 



0.4 0.6 
orbital phase 
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Fig. 1. Mean RVs obtained by direct measurement. Donor ve- 
locities were obtained averaging the Doppler shifts for the 
MgII4481A, FeI4957A, Sill 6347 A, and Sill 6371 A lines. 
Hel lines were used to calculate gainer velocities. 



o 



1.0 



0.5 



0.0 



1 1 1 1 1 1 


1 1 

Composite Spectrum . 

Donor 




Gainer 


i i i i i 



4000 4100 



4200 4300 4400 
Wavelength [A] 



4500 4600 



Fig. 2. Separated spectra of DQ Vel. The composite spectrum 
(top) at (f> = 0.81 10 is used here as a reference from the list in 
Table 1 . Final S/N for the decomposed spectra are 97 (middle) 
and 267 (bottom) for the donor and gainer star. The spectra were 
offset for clarity. 



using a narrow CCF (Hilditch 200 li, we applied separating and 
CCR processes only to the CORALIE data. The main reason 
for this is that lines on the template spectrum and the observed 
spectra will have similar shapes and widths, which decreases the 
broadening of the CCF. 

Using previously calculated CORALIE RVs, we shifted all 
spectra to the donor rest-frame and combined these to create a 
first-guess spectral template for the donor star. This template 
is an orbital-cycle averaged single-line spectrum. This proce- 
dure attenuates the spectral features of the remaining component 
(gainer). Subtracting this donor template from all spectra and 
shifting the resulting spectra to the gainer rest-frame, we com- 
bined the data again to obtain an orbital-cycle-averaged template 
for the gainer star. We repeated this procedure until the residuals 
of one component contribution disappeared from the template 



Table 3. Calculated RVs for DQ Velorum using ID CCR. RV d 
and RV S are the mean RV for the donor and gainer star. 



BJD-2450000 



RV d (kms- 1 ) /fVgOons" 1 ) 



4612 


5771 





0247 


62 


28 


± 


0.64 


-58 


99 


±0.63 


4612 


6069 





0296 


66 


37 


+ 


0.59 


-6 


90 


±0.49 


4612 


6367 





0345 


81 


65 


± 


0.46 


2 


12 


±0.32 


4819 


7745 





0843 


139 


21 


± 


0.60 


-2 


74 


±0.42 


4564 


6168 





1408 


157 


54 


± 


0.69 


-57 


80 


±2.48 


4564 


6466 





1457 


163 


68 


± 


0.53 


-97 


07 


±0.42 


4564 


6764 





1506 


164 


80 


± 


0.52 


-78 


05 


±0.30 


4936 


6295 





2932 


161 


90 


± 


0.42 


-68 


49 


±0.41 


4936 


6594 





2981 


160 


04 


± 


0.61 


-57 


92 


±0.42 


4565 


5829 





2996 


150 


65 


± 


0.43 


-67 


78 


± 1.03 


4936 


6892 





3030 


153 


33 


+ 


0.62 


-56 


34 


±0.53 


4565 


6127 





3045 


143 


18 


+ 


0.52 


-46 


62 


±0.54 


4565 


6425 





3094 


146 


92 


+ 


0.54 


-41 


47 


±0.89 


4937 


6087 





4542 


3 


89 


+ 


0.54 


-52 


49 


±0.42 


4937 


6386 





4591 


-17 


44 


+ 


0.77 


9 


33 


±0.57 


4937 


6684 





4640 


-18 


72 


+ 


0.53 


-21 


33 


±0.37 


4609 


5451 





5263 


-95 


82 


+ 


0.50 


-23 


60 


±0.73 


4609 


5750 





5312 


-99 


43 


+ 


0.52 


26 


70 


±0.52 


4968 


5220 





5358 


-95 


47 


+ 


0.54 


18 


83 


±0.48 


4609 


6258 





5395 


-108 


24 


+ 


0.51 


34 


00 


±0.44 


4968 


5507 





5405 


-97 


11 


± 


0.57 


-44 


20 


±0.33 


4822 


7697 





5767 


-131 


58 


+ 


0.58 


34 


14 


±0.37 


4938 


6259 





6214 


-157 


34 


+ 


0.52 


29 


28 


±0.45 


4938 


6701 





6287 


-160 


61 


+ 


0.49 


15 


12 


±0.27 


4938 


6999 





6336 


-159 


57 


+ 


0.49 





25 


±0.24 


4610 


5343 





6889 


-184 


89 


+ 


0.46 


37 


11 


±0.59 


4610 


5641 





6938 


-183 


61 


+ 


0.51 


31 


42 


±0.43 


4610 


5940 





6987 


-183 


13 


+ 


0.50 


54 


30 


±0.38 


4969 


5085 





6980 


-180 


27 


+ 


0.64 


75 


60 


±0.37 


4969 


5378 





7028 


-177 


57 


+ 


0.50 


49 


58 


±0.34 


4562 


5807 





8061 


-141 


67 


+ 


0.44 


56 


73 


±0.42 


4562 


6105 





8110 


-146 


29 


+ 


0.49 


58 


45 


±0.38 


4562 


6403 





8159 


-144 


55 


+ 


0.47 


53 


10 


±0.37 


4611 


5221 





8512 


-112 


38 


+ 


0.57 


55 


44 


±0.39 


4970 


4967 





8604 


-96 


28 


+ 


0.41 


3 


68 


±0.49 


4970 


5253 





8651 


-84 


37 


+ 


0.59 


55 


73 


±0.37 


4611 


6078 





8653 


-96 


31 


+ 


0.61 


44 


25 


±0.49 


4611 


6377 





8702 


-97 


31 


+ 


0.59 


27 


13 


±0.37 



spectra of the other. According to Gonzalez & Levato (2006i, 
no more than 5-7 iterations are necessary considering that with 
each iteration the residuals are reduced by a factor 1 In, where n 
is the number of observed spectra. Figure 2 shows a portion of 
our separated spectra after six iterations. 

We applied a ID Fourier CCR process to the observed spectra 
and previously obtained templates using the IRAF fxcor task. 
For binaries this task works as follows: we select the template of 
one component (which is by construction at rest in the barycen- 
tric frame), and convolve it with each of the observed spectra in 
which the spectral lines move according to the relative motion 
of the binary. The maximum of the CCF will be, in this case, 
the shifted velocity of the spectrum relative to the template. The 
process is then repeated using the other template. In this way, we 
obtain the relative velocities of each spectrum with respect to the 
donor and gainer templates. 

During the CCR process and to finally remove the weak modula- 
tion found on the CORALIE spectra, we applied a Hanning-type 
filter to the template and the observed spectra after the data were 
transformed and prior to correlation. This Hanning-filter func- 
tion is used to attenuate the Fourier components over a specific 
wavenumber range to remove remaining high-frequency noise 
that was not removed by the continuum subtraction. 
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Table 4. RVs for DuPont/echelle and MPI/FEROS spectra not 
included in the CCR process. 



BJD-2450000 


<t>o 


tfVdCkms- 1 ) 


tfV g (km s- 1 ) 


4631.4735 


0.1309 


139.40 


+ 


1.33 


-19.78 


±8.38 


4631.4854 


0.1329 


143.62 


+ 


1.37 


-53.00 ± 


11.23 


4966.4639 


0.1975 


177.57 


± 


0.89 


-66.54 


±4.28 


4966.5077 


0.2047 


185.23 


± 


1.72 


-36.84 


±7.14 


5283.5531 


0.3214 


146.55 


+ 


3.13 


-69.58 


±4.41 


4967.5020 


0.3681 


123.41 


+ 


1.00 


-27.20 ± 


10.91 


5285.5668 


0.6524 


-169.66 


+ 


1.79 


63.94 


±4.26 



200 F 




-200 t 



0.0 



0.2 



0.4 0.6 

orbital phase 



0.8 



.0 



Fig. 3. RV curves and best sinus fit for DQ Vel components. Each 
RV was measured from Table 3 and Table 4 taking the mean 
value in orbital phase bins of 0.05. 

Table 5. Circular-orbit solutions obtained for the best sinus fits. 



Parameter 


Donor star 


Gainer star 


K(km s- 1 ) 


185.436 ±4.963 


-57.556 ±5.536 


y (km s -1 ) 


1.858 ±3.712 


-8.745 ± 4.066 


r.m.s (km s _1 ) 


14.80 


14.07 



Notes. K is the semi-amplitude of RV curves and y is the systemic ve- 
locity. Data included here were obtained from the CCR process and 
direct measurements. 



From the donor template, we calculated RVs independently in 
two spectral ranges at 4500-4800 A and 4900-5600 A, where we 
found more spectral lines from the donor star. A third measure- 
ment was made in both spectral ranges simultaneously. Final 
RVs were obtained by taking the mean between these three mea- 
surements. We applied the same procedure using the gainer tem- 
plate in the 4350-4750 A and 4900-5300 A ranges. Table 3 shows 
our final RVs determined from the CCR process and their corre- 
sponding r.m.s. errors. 

In Table 4 we show the RVs for the remaining spectra that were 
not included in the CCR process. 

To use all available RV measurements, we decided to obtain the 
mean RV value in orbital phase bins of 0.05. We fitted the RV 
curves with a sinus function. Best fits and circular-orbit solutions 
are shown in Figure 3 and Table 5. The different y obtained for 
the gainer suggests that their lines are influenced by circumstel- 
lar material. 



Based on these results, we determined the projected semi- 
major axes and minimum masses for donor and gainer star using 
the expressions by Hilditch (1200 11 1 for a circular orbit: 

a d ,g sin i = (1.9758 x 10" 2 )^ d „P o [R Q ] (2) 

a d sin i = 22.288 + 0.596 [R B ] (3) 

fl g sin / = 6.917 + 0.665 [R Q ] (4) 

m d , g sin 3 i = (1.0361 x 10~ 7 )(K d + K g ) 2 K g4 P D [M Q ] (5) 

m d sin 3 i = 2.141 + 0.244 [M ] (6) 

m g sin 3 i = 6.901 ± 0.460 [M ] , (7) 

where a d and a g are the donor and gainer semi-major axes mea- 
sured in solar radii (R Q ), m d and m g are the donor and gainer 
masses measured in solar masses (M ), K& and K g are the semi- 
amplitudes of the RV curves expressed in kms -1 , P a is the or- 
bital period in days, and i the orbital inclination of the system. 
It can be shown that independently of the orbital inclination, 
equation 5 yields 



m d _K g _ 
m g Ka 



0.31+0.03. 



(8) 



which represents the spectroscopic mass ratio q of the system. 
Note that for this q value we assumed that the observed RV vari- 
ations of the gainer's lines represent the motion of the gainer. 
If we consider that the donor star fills its Roche lobe, we can ap- 
proximate the donor radius (Rd) to the size of the donor Roche 
lobe according to Eggleton (119831 1: 



r L,d 



a 



QA9q 



2/3 



0.69^ 3 + ln(l +qV 3 ) 



where a is the semimajor axis of the relative orbit. 
For q = 0.31 + 0.03, we obtain 

R d = (0.269 + 0.019) a. 



(9) 



(10) 



For an eclipsing binary with an orbital inclination 70° < i < 90° 
we can estimate a donor radius between (7.856 + 0.604) R Q < 
Rd ^ (8.366 ±0.643) R G , which indicates an oversized star based 
on the previous mass value. 

We can also determine the maximum value for the equatorial 
rotational velocity of the donor assuming it is rotating syn- 
chronously: 



2nR 



V r ,d 



Po 



d 69.606 + 5.349 km s" 1 



(ID 



for which we used the maximum donor radius obtained above. 

3.2. Spectral type classification and analysis of spectral 
features 

The donor star template shows several absorption metal lines in- 
cluding Fel - II, Mgl - II, Till, CrII, and Sill. Some lines like 
Fell 4233.2 A, Fel 4383.5 A, Mgl 5 167.3 A, Mgll 5528.4 A, and 
Sill 5688.8 A were not clearly identified in the original spectra. 
We compared our donor template with a grid of synthetic spec- 
tra in the region 5200-5350 A which is deprived of hydrogen and 
helium lines but has several metal lines. To determine the grid of 
synthetic fluxes we used atmospheric models computed with the 
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Fig. 4. Calibration curves for measured fwhm (A) of 
Mgl 4703.0 A, Mgll 448 1.3 A, Till 4780.0/4805.1 A and 
Sill 637 1.3 A lines plotted against v r sin 



line-blanketed LTE ATLAS9 code (Kurucz IT9931) . which treats 
line opacity with opacity distribution functions (ODFs). The 
Kurucz models are constructed assuming plane-parallel geome- 
try and hydrostatic and radiative equilibrium of the gas. The syn- 
thetic spectra were computed with the SYNTHE code (Kurucz 
1993). Both codes, ATLAS9 a nd SY NTHE, were ported un- 
der GNU Linux by Sbordone (2005 ) and are available online 
5 . The atomic data were taken from Castelli & Hubrig 120041 
\ . The theoretical models were obtained for effective tempera- 
tures from 6000 to 10000 K with steps of 100 K and for sur- 
face gravities from 2.0 to 4.5 dex with steps of 0.1 dex. Solar 
and 0.5 dex higher metallicities were taken into account. The 
grid of synthetic spectra was calculated for five different rota- 
tion velocities, v r sinz' = 0, 25, 50, 75 and 100 kms -1 . We sub- 
tracted the donor template from every grid spectrum and anal- 
ysed the residuals of the resulting spectra. The best model is 
obtained for an A-type star with 7^ = 9400 K, loggd = 3.1, 
v r sin/ = 75km s _1 and solar metallicity. To obtain a more re- 
stricted value for the rotational velocity, we decided to inter- 
polate between the different values of rotation velocities gener- 
ated in the grid of synthetic spectra. We selected four synthetic 
spectra from the grid with T = 9400 K, logg = 3.1 and ve- 
locities of 25, 50, 75, and 100 kms _I . We measured the full 
width at half maximum (fwhm) for the isolated and unblended 
metal lines Mgl 4703 .0 A, Mgll 448 1 .3 A, Till 4780.0/4805 . 1 A, 
and Sill 637 1.3 A. We plot in Figure 4 the fwhm-Vj sin i cali- 
bration curves. We calculated the fwhm for the same spectral 
lines in the donor template and determined a mean velocity of 
v r ,d sin i — 65 + 4 km s _1 , which agrees with our previous result 
of equation 1 1 . 

A first inspection of spectral features in the gainer template sug- 
gests a B-type star with marked Balmer and helium lines like 
FM4026/ 4471/4921/5015/5875/ 6678Aand also a few and 
weak metal lines such as Fell 4303 / 5 169/ 5276 A, Mgll 448 1 A, 
and SHI4133 A. We also observed asymmetric profiles in the he- 
lium and Balmer lines after the separation, as shown in Figure 7. 
In the current DPV scenario, the donor star fills its Roche lobe 
and transfers mass to a rapidly rotating gainer star via gas stream. 




3 wwwuser.oat.ts.astro.it/atmos/ 

4 http://wwwuser.oat.ts.astro.it/castelli/grids.html 



Fig. 5. Donor light contribution with respect to the total light cal- 
culated from the light curve model. Each curve shows the donor 
contribution for a specific spectral range around a spectral line. 
The selected spectral lines (and ranges) are from top to the bot- 
tom: Hel 6678 A, Ha, HeI5875 A, H/3, HeI4471 A, and Hy. LI 
and L2 are the stellar fluxes and Ld is the disc contribution. 



Therefore, the analysis of the gainer spectra requieres special 
attention to check for the possible presence of circumstellar 
material. To study the spectral-orbital behaviour of the gainer, 
we decided to remove the donor contribution from the origi- 
nal spectra using the previously found donor template. To de- 
rive the contribution of the donor to the total light along the or- 
bital phase, we modeled the light curve in the V-band with the 
previously calculated donor temperature and mass ratio and as- 
sumed that the donor light adds to the total light. Under these 
assumptions, the total light includes two stellar components (LI 
and L2) plus a circumstellar disc component (Ld). More details 
about the photometric model are explained in section 3.4. We 
calculated the donor factor contribution to the total light for spe- 
cific spectral lines. These factors were determined taking into 
account the orbital changes in the proyected area of the donor 
as well as the wavelength dependence on the different colour 
terms of donor and gainer (see Fig. 5). From each observed spec- 
trum, we subtracted the donor template multiplied by the corre- 
sponding factor for the specific wavelength range. The resulting 
spectra were normalised to the new continuum level. Selected 
donor-subtracted Ha and H/5 profiles are shown in Figure 6. The 
two profiles are composed of a central absorption surrounded 
by weak blue/red asymmetric wings, whose velocities reach up 
to about ~ 350 kms 1 and are phased with the gainer's RVs. 
Similar behaviour was found in Hy profiles. We detected a weak- 
ness of central absorption during the orbital phases <p - 0. 1 -0.2 
and <p = 0.6 - 0.7 also observed in the helium lines and proba- 
bly associated to a filling emission from a higher local tempera- 
ture interaction regions such as hot or bright spots, visible during 
those orbital phases. 

To gain some insight into the gainer temperature we calcu- 
lated equivalent widths (hereafter E W) for the Hel 447 1 A and 
Mgll 4481 A lines for all donor-subtracted spectra. Mennickent 
at al. d2012bl l determined the EWs ratios R = EW44ji/EW44%\ 
for B-type stars using a grid of synthetic spectra (log g = 4.0, 
Vturb - 2km s _1 and solar metallicity) and plotted it against 
temperatures. These results were used to determine the tem- 
perature of the gainer. Similarly, we calculated R in all donor- 
subtracted spectra and used the mean value to estimate a tem- 
perature T g ~ (15 000 + 1000) K for the gainer star. We show the 
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Fig. 6. Selected donor-subtracted Ha (left) and H/3 (right) pro- 
files along the orbital period. Orbital phases are shown on the 
right side of the left plot. Same values were used in the right 
panel. 



variability of £W447i and along the orbital phase in Figure 8. 



3.3. Orbital and long light curves 

We separated multiperiodic V, I, J, and K light curves with a 
Fourier decomposition technique (for more details see Demircan 
1 19791 Mennickent et al. 120 1 2al > and found two main Fourier 
components of frequencies /i and /i, the sum of which (includ- 
ing their most important harmonics) is the best representation 
of the light curves, /i is higher and associated with a period of 
6.083299 d, which represents the orbital variability of the binary 
(see Fig. 9). We find the following orbital cycle ephemeris for the 
light curve minimum: 



HJD min , orh = 2453407.60(2) + 6.083299(7) x E , 



(12) 



which agrees with the ephemeris from van Houten ( 1950 1. 
The frequency is related to a period of 188.7 d and is rep- 
resented by a smooth variability whose amplitude in V-band 
is ~ 2% with respect to the total luminosity (see Fig. 10). The 
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Fig. 7. Asymmetric helium and Balmer profiles found in the 
gainer template after separating the spectra. 
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Fig. 8. EW4471 (bottom) and R = W4471/W4481 (top) plotted 
against orbital phase. 



ephemeris for the maximum of the long-term variability is best 
represent by 



HJD 



max, long 



2453437.2(16) + 188.7(2) x E . 



(13) 



A visual inspection of the orbital light curves shows an eclipsing 
variable with well-defined eclipses and noticeable proximity 
effects. Clearly, at least one of the components is tidally dis- 
torted. Differences in the depths of the minima suggest different 
surfaces temperatures. The deepest primary minimum on the 
V-band indicates that the bluer, hotter, and smaller gainer star 
is being eclipsed. Its light contribution is even less dominant 
on the infrared bands, supporting our previous results from the 
spectroscopic studies. 



3.4. Light curve modelling 

To derive the main physical parameters for the stellar compo- 
nents we fitted the V-band light curve with the inverse-problem- 
solving method based on the simplex algorithm. Main binary 
elements used in the model and the light-curve-fitting procedure 
are described in Djurasevic (1992] [1996). We assumed a semi- 
detached configuration with the donor filling its Roche lobe. 
After unsuccessful attempts to model the system without a disc, 
we finally adopted a configuration that included an optically 
thick accretion disc around the gainer star. The disc is charac- 
terised by its radially dependent temperature T^dr), its radius 
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Fig. 9. Separated VIJK orbital light curves during one season in 
2008 phased with the orbital period. 
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Fig. 10. DQ Vel ASAS V-band light curve after separating the 
short-term (top) and long-term (bottom) variability. 



iodise, its thickness at edge d e and at centre d c , and by assuming 
that it is in physical and thermal contact with the gainer. The 
temperature distribution of the disc is given by 



T asd r) - Tdiscc + (To- 7discJ 



1 



(14) 



where 7di SCe is the temperature of the disc at the edge, a 1 is 
the temperature exponent, and (T g ,R g ) represent the tempera- 
ture and radius of the gainer star. 

The model also includes two active regions, hot spot (hs) and 



bright spot (bs), with higher local temperatures, which are lo- 
cated at the edge of the disc. Including these regions signifi- 
cantly improved the fit. They are represented by their tempera- 
tures (7h s ,bs), angular dimensions (Shs.bsX an d longitudes (/lhs,t>s), 
which are free parameters in the model. The longitudes /lh S ,bs are 
measured clockwise in the range - 360 degrees and are viewed 
from the direction of the +z-axis, which is orthogonal to the or- 
bital plane. More details about the disc model are described in 
Mennickent et al. l2012al 

To restrict the number of free parameteres we fixed the donor 
temperature Td = 9400 K and mass ratio q = 0.31, whose val- 
ues were reliably derived from the previous spectroscopic anal- 
ysis (Sect.3.1 and 3.2). We also fixed the gravity-darkening and 
albedo coefficients of the components to the values f3d, g = 0.25 
and Ad, g = 1.0 according to the von Zeipel law for radiative 
shells (von Zeipel 1924). Considering that the donor star fills its 
Roche lobe, we assumed this star to rotate synchronously with 
a rotation coefficient fd = 1 .0. We modeled the ASAS V-band 
light curve with the gainer star in synchronous-rotation regime 
(/ g = 1.0) and also in a critical-rotation scenario with a non- 
synchronous rotation coefficient f g = 12.8. We observed no sig- 
nificant difference in the two fitting results, revealing a minor 
effect of the gainer rotational velocity on the system parameters. 
The V-band light curve fitting with the gainer in synchronous 
rotation is shown in Figure 1 1 . Table 6 shows the final parame- 
ters of the fitting. The main physical properties of the gainer star 
such as mass (M g = 7.3 + 0.3 M Q ), radius (3.6 ± 0.2 R Q ), and 
gravity (log g g = 4.2+0.1) agree with our previous spectroscopic 
result for the spectral classification (B-type star). However, we 
find a temperature value for the gainer T g = 18 500 + 500 K 
different from the previously obtained T g = 15 000 + 1000K 
using the equivalent width ratio between the spectral features 
R = EWwiil EWaazx- This difference is probably associated with 
the disc contribution to the spectra. 

As was mentioned before, the best-fit model for DQ Vel includes 
an optically thick disc around the gainer star. From Figure 1 1, 
we estimate a disc contribution of about 7% of the total light at 
quadrature phases. It has a concave shape with higher thickness 
at its edge than at its centre and is extended. It fills (89 + 3)% 
of the critical gainer Roche lobe radius with a temperature of 
6580 ± 300 K at the outer radius, as shown in Table 6. According 
to the model, the hot spot -where the material coming from the 
donor impacts the disc- is located at longitude 329.1° + 7.0° 
and has a temperature of 7h s = 9500 + 750K, which is close 
to the donor temperature found previously. This is probably the 
reason why we do not observe emission from the hot spot on 
the composite spectra. Because the disc is extended, the distance 
between LI and the outer radius of the disc is short and not suffi- 
cient to produce a substantial acceleration of the gas stream and 
thus reach a higher temperature at this region. 

3.5. Possible pulsations of the B-type gainer star 

As suggested in Sect. 3.1, the scatter in the gainer RVs might be 
associated with pulsations of the B-type gainer star. The main 
reason to consider this scenario is the location of the gainer in 
the H-R diagram, as shown in Figure 12. According to its mass 
and luminosity, the gainer is located well within the instability 
region of slowly pulsating B-stars (SPBs). These pulsators are 
typically B2 to B9 stars of intermediate mass between 3 to 7 
M . SPBs show multiperiodic light variations of the order of 
days produced by high-order g-modes. The typical photometric 
amplitude variations are lower than 0. 1 magnitude while the ra- 
dial velocity variations can reach up to 15 km s~' (De Cat 2002 
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Table 6. Analysis results of the DQ Vel V-band light-curve. 



Quantity Remarks 



n 


583 


Number of observations. 


1.(0 - C f 


0.3984 


Final sum of the residuals square between observed (LCO) and synthetic (LCC) light-curves. 




0.0258 


Root-mean-square of the residuals. 


f'° 


82.5 ±0.2 


Orbit inclination (in arc degrees). 


Fdisc = ^disc/^yc 


0.89 ± 0.03 


Disc dimension factor (the ratio of the disc radius to the critical Roche lobe radius along the y-axis) 


^disc[K] 


6580 ± 300 


Disc-edge temperature. 


rfc[«orb] 


0.019 ± 0.005 


Disc thicknesses at the edge in the units of the distance between the components. 


d c [a mb ] 


0.011 ±0.005 


Disc thicknesses at the centre in the units of the distance between the components. 




4.5 ±0.3 


Disc temperature distribution coefficient. 


F g 


0.271 ±0.005 


Filling factor for the critical Roche lobe of the gainer in synchronous rotation regime. 


r g [K] 


18 500 ±500 


Temperature of the gainer. 


Ahs = ^hs/Tdisc 


1.45 ±0.1 


Hot-spot temperature coefficient. 


0hs[°] 


16.0 ±2.0 


Hot-spot angular dimension. 


4 S [°] 


329.1 ±7.0 


Hot-spot longitude (in arc degrees). 


0rad[°] 


-19.7 ±5.0 


Angle between the line perpendicular to the disc edge and the direction of the hot spot. 


j4bs = T bs /T iisc 


1.39 ±0.1 


Bright-spot temperature coefficient. 




50.1 ±5.0 


Bright-spot angular dimension. 




142.7 ± 9.0 


Bright-spot longitude (in arc degrees). 


n g 


8.45 ± 0.02 


Dimensionless surface potentials of the gainer. 


fid 


2.49 ± 0.02 


Dimensionless surface potentials of the donor. 


M e [M ] 


7.3 ±0.3 


Mass of the gainer (in solar masses). 


M d [M ] 


2.2 ±0.2 


Mass of the donor (in solar masses). 


/? g [i? ] 


3.6 ±0.2 


Mean radius of the gainer (in solar radii). 


#d [/?o] 


8.4 ± 0.2 


Mean radius of the donor (in solar radii). 


!°ggg 


4.2 ±0.1 


Logarithm of the gainer's effective gravity. 


loggd 


2.9 ±0.1 


Logarithm of the donor's effective gravity. 


ML 


-3.1 ±0.2 


Absolute bolometric magnitude of the gainer. 


<> 


-1.9 ±0.1 


Absolute bolometric magnitude of the donor. 




29.7 ±0.3 


Orbital semi-major axis. 


^disc L^o] 


12.9 ±0.3 


Disc radius (in solar radii). 




0.6 ±0.1 


Disc thicknesses at the edge in solar units. 


d c [R e ] 


0.3 ±0.1 


Disc thicknesses at the centre in solar units. 



Notes. Fixed parameters: q = M d /M g = 0.31 -mass ratio of the components, T A = 9400K -temperature of the less-massive donor, F d = 1.0 
-filling factor for the critical Roche lobe of the donor, / g _d = 1.00 -synchronous rotation coefficients of the system components, yS g _d = 0.25 
-gravity-darkening coefficients of the components, A g d = 1.0 -albedo coefficients of the components. 



De Cat & Aerts l2002l l. 

To find possible stellar pulsations, we carried out a Fourier fre- 
quency analysis on the gainer RVs and also on the O-C residuals 
of the V-band light curve after subtracting the orbital and long 
variabilities (see second panel in Figure 11). 
The Fourier frequency analysis was performed with the software 
Period04E which is based on a discrete Fourier transform algo- 
rithm and is applied to the spectral analysis of unevenly sampled 
data. 

We computed the Fourier transform (hereafter FT) for the gainer 
RVs data set. We found a significant peak of 0.16 d _1 , which 
corresponds to the binary orbital frequency. In addition, several 
marginal peaks appear at higher frequencies, but the sampling of 
our scarce data set is not sufficient to unambiguously determine 
any periodicities above the noise level. 

Using the larger sample of 516 data points obtained after sub- 
tracting the orbital and long variability from the V-band light 
curve, we calculated the FT for this data set. In the top and bot- 
tom panels of Figure 13 we show the FT for the whole sample 
and its corresponding spectral window up to the Nyquist fre- 
quency (//v ~ 0.25 d" 1 ). To distinguish the significance of a 
frequency peak we used a 4cr„ , je criterion, which resulted in a 
99.9% confidence limit (Kuschnig et al. l 19971 ). To determine the 
Cnoise in the Fourier spectrum we calculated the noise level by 
averaging the amplitudes around each frequency along the whole 

5 http://www.univie.ac.at/tops/Period04/ 



frequency domain. As shown in Figure 13, the solid white line 
in the FT corresponds to this mean noise level. The dotted line 
indicates a significance threshold of four times the cr no j se . After 
determining the highest frequency peak above 4cr no(je in the FT 
of the original data set, we prewhitened the data, removing the 
sinusoid corresponding to this frequency. Then, we calculated 
the FT of the residuals and found the next highest peak. We re- 
peated this procedure until no new peaks with significance above 
4cr„ OIje could be detected. In Table 7 we list the seven detected 
frequency signals. According to their S/N, only two frequencies 
if i and fx) are considered significant. In the middle upper panel 
of Figure 13 we show the FT of the residuals after removing f\ . 
Their corresponding cr no!je and 4<x„ , J( , levels are shown. In the 
middle lower panel the prewhitened FT after removing f\ and fx 
is displayed. No new significant peaks appear in the power spec- 
trum after removing f\ and fx. The frequencies f\ = 0.005 d _1 
and fa — 0.166 d _1 correspond to the long and orbital photomet- 
ric variabilities, which are still present in the data. The origin of 
the significant frequency fx = 0.194d~' corresponding to a pe- 
riod of 5. 14 d is unknown and could be interpreted as a pulsation 
of the SPB type. 

We independently computed a simple Monte Carlo simulation 
to check the significance of fx. To do this, we generated 1000 
samples of synthetic random data with a Gaussian distribution. 
For each sample, we generated an amount of points equal to the 
original number of 516 O-C data points and used the same time 
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Fig. 11. Panels from top to bottom: Observed (LCO) and syn- 
thetic (LCC) light-curves of DQ Vel obtained by analysing pho- 
tometric observations; final O-C residuals between the observed 
and optimum synthetic light curves; fluxes of donor, gainer and 
of the accretion disc, normalised to the donor flux at phase 0.25; 
the views of the optimal model at orbital phases 0.25, 0.50 and 
0.75, obtained with parameters estimated by the light curve anal- 
ysis. 




RR Lyr 




Fig. 12. Location of DQ Vel including the errors of stellar com- 
ponents on the H-R diagram. Some stellar pulsation instability 
regions are indicated. The evolutionary tracks were taken from 
Bressan et al. [T993l for z=0.02. 

Table 7. Significant and marginal frequencies detected in the 
Fourier spectrum of the (O-C) DQ Vel photometric data. 



ID 


frequency 


amplitude 


phase 


S/N 




(d- 1 ) 


(mma) 






h 


0.005 


24 


0.65 


6.04 


h 


0.194 


17 


0.43 


5.18 


n 


0.166 


11 


0.82 


3.67 


n 


0.086 


11 


0.38 


3.67 




0.030 


11 


0.37 


3.56 


n 


0.148 


9 


0.14 


3.08 


fi 


0.102 


9 


0.74 


2.95 



Notes. 'According to their S/N, these frequencies are considered 
marginal detections. 

sampling. Then, we calculated the FT for each data set in the 
same frequency range as in the previous analysis, and found the 
amplitude of the highest frequency peak. In 99.9% of cases we 
do not find any amplitude peak higher than 17 mma. This is an 
additional independent test indicating the significance of the fj 
detection. 

3.6. Distance to DQ Vel 

It is well known that eclipsing double-lined binaries can be used 
as standard candles to calculate accurate distances for Local 
Group galaxies (Kang et al. 2007, Clausen 2004). The most of- 
ten used method for Milky Way binaries allows one to derive 
the distance modulus for either component of a binary (Clausen 
120041 1. In the V-band we have 



(m d , g - M djg ) = 5 log(R d JR e ) + (m d , g - A v ) 

- M bol0 + lOlog(r d , g /r ) + BC d o 



(15) 



The absolute radii of stars and effective temperatures together 
with the visual apparent magnitudes of the components (m d , m g ), 
interstellar absorption (A v ), and bolometric corrections (BC) 
are needed. Considering the apparent magnitudes on the light 
curves including the disc contribution, we determined the ap- 
parent magnitudes for the donor and gainer star using the V- 
band light curve and the theoretical fluxes from the photomet- 
ric model (Fig. 11). From this last we measured the total flux 
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Fig. 13. Top panel: FT of the 516 data points obtained from the 
(O-C) DQ Vel photometric data. Middle upper panel: FT of the 
residuals after prewhitening the original data by the first signif- 
icant frequency (/i). Middle lower panel: FT of the residuals 
prewhitened by the two significant frequencies f\ and listed 
in Table 7. The white solid line corresponds to cr„ 0!je and the dot- 
ted line corresponds to 4cr noiJ( ,. Bottom panel: Spectral window 
for the whole data set. 



of the system (/ t ) and the individual contributions (fa, f g , /di sc X 
where f t = fa + f g + /disc in the cuadrature phases (p = 0.25 and 
also (f> = 0.75. The total flux at these orbital phases is related 
to the highest observed apparent magnitude of the system in the 
light curve m t (V) = 10.977. Thus, the apparent magnitude of the 
donor can be obtained from 



m,\ - f'<h = -2.5 log) ^ 



11.578 + 0.01 



(16) 



The same procedure was applied to the gainer star to obtain 
m g = 12.097 + 0.01. The bolometric corrections were obtained 
from the calibrations by Flower (1996 1 using log T g = 4.267; 
BC g = (-1.76 + 0.16) for the gainer and log T d =°3.973 and 
BCd = (-0.13 + 0.05) for the donor. We determined the inter- 
stellar absorption A v = 1 .02 + 0.02 using the algorithm from 



Table 8. Comparison of the main stellar and disc parameters for 
DQ Vel and V393 Sco. 



Parameter 


DQ Vel 


V393 Sco 


Stellar Components 


r g [K] 


18 500 ±500 


16 600 + 500 


r d [K] 


9400 ± 100 


8600 ± 600 


in 


82.5 ± 0.2 


80.0 ± 0.2 


M, [M ] 


7.3+0.3 


7.8 + 0.2 


M d [M Q ] 


2.2 + 0.2 


2.0 ± 0.2 


Rg [«0] 


3.6 ± 0.2 


4.4 ± 0.2 


Ri [R e ] 


8.4 + 0.2 


9.4 ± 0.3 




4.2 + 0.1 


4.0 + 0.1 


loggd 


2.9 + 0.1 


2.8 + 0.1 


"orb [^ol 


29.7 ± 0.3 


35.1pm0.3 


Disc 




0.89 ± 0.03 


0.55 ± 0.04 


^disc L^ol 


12.9 + 0.3 


9.7 ± 0.3 


ck [R Q ] 


0.6 + 0.1 


1.3 + 0.3 


d c [R ] 


0.3+0.1 


2.1+0.4 



Notes. Results for V393 Sco were taken from Model A in Mennickent 
et al. d2012ab . Parameters for DQ Vel are taken from this work (Table 
6). 

Cardelli et al. (1 19891 1 that is available online^. 
Equation 15 yields 



(m g - M g ) = 12.46 + 0.23 
(m d - M d ) = 12.47 ± 0.09 . 



(17) 
(18) 



Finally, we estimated the average distance from donor and gainer 
distance modulus: 



distance(pc) = io« m - M >° +5 >/ 5 
distance(pc) — 3112+175. 



3. 7. Comparison with V393 Scorpii 



(19) 



During recent years a comprehensive study of the DPV sys- 
tem V393 Sco has been carried out by Mennickent at al. (2010 
I2012al 120 1 2bt to study the nature and evolution of the system. 
The studies include photometry (VIJK bandpasses and V-ASAS 
light curves) and high-resolution spectroscopy in the optical, 
near-infrared, and ultraviolet wavelength ranges. V393 Sco is a 
galactic (d = 523 pc) bright DPV composed of an A-type donor 
(r d = 8600 ± 600 K) and a B-type gainer (T g = 16600 ± 500 K) 
surrounded by a massive disc. The system shows a complex vari- 
ability of spectral features along the orbital and long cycles in- 
cluding intense double-peaked emission in optical Balmer lines. 

We noted strong similarities of the physical properties of V393 
Sco and DQ Vel, as can be seen in Table 8. Physical parameters 
of the stellar components such as masses, radii, temperatures, 
gravities, and also orbital parameters (binary separation and in- 
clination) are similar. However, a detailed comparison showed 
significant differences in the optical spectral features of the two 
systems. We suggest that these differences are associated with 
the physical properties of the accretion discs and thus with the 
mass transfer rates. As shown in Table 8, the disc in DQ Vel 
seems to be more extended and colder than that of V393 Sco. 
Moreover, the hot and bright spots in V393 Sco are hotter than 
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in DQ Vel. From the geometric point of view, DQ Vel's disc is 
concave and thick at the outer border, in contrast to the mod- 
elled disc in V393 Sco, which is convex and probably has most 
of its mass concentrated close to the gainer. We suggest that the 
main differences of the discs are related to different evolution- 
ary stages. The high mass-transfer rate of V393 Sco would form 
a massive accretion disc and lead the gainer to rotate critically. 
DQ Vel instead could be an older system in a state of lower mass 
transfer where the gainer had time to slow down, allowing the 
formation of an extended disc. A detailed analysis of the evolu- 
tionary stage of DQ Vel remains to be conducted in the future. 



4. Summary 

Our analysis of VIJK (REM-telescope) and V (ASAS) photo- 
metric data togheter with the study of 46 optical high-resolution 
echelle spectra along several orbital cycles allowed us to esti- 
mate the main physical properties for the components of the 
DPV system DQ Vel. 

We separated the composite spectra using the Doppler shifts of 
known spectral lines to obtain two single-lined templates, which 
were used to derive the RV curves through a Fourier cross- 
correlation process, and obtained the spectroscopic mass ratio 
(q = 0.31+ 0.03) of the system. A comparison of the templates 
with a grid of synthetic spectra enabled us to estimate the effec- 
tive temperature = 9400 + 100 K for the donor star. We cal- 
culated the fractional light contribution of the donor star to sub- 
tract it from the spectra and thus analyse the gainer spectral fea- 
tures. The donor-subtracted spectra along the orbital cycle show 
a variability of the Balmer and helium profiles, suggesting the 
presence of circumstellar material. The spectroscopic mass ra- 
tio and donor temperature were used to carry out a V-band light 
curve fitting and to obtain the main binary elements such as the 
masses M g = 7.3 + 0.3 M and M d = 2.2 ± 0.2 M , and the radii 
R g = 3.6 + 0.2 Rq and/?d = 8.4 + 0.2 R for gainer and donor star 
together with the orbital separation (a or b = 29.7+0.3 R e ) and the 
inclination angle (;'° = 82.5 ± 0.2) of the system. The best model 
suggest a semi-detached system comprising an A1III donor star 
(r d = 9400+100 K) and a B3V gainer star (T g = 18 500+500 K) 
with an extended accretion disc around the gainer. The model 
includes two active regions (hot and bright spots) located at the 
outer edge of the disc. The temperatures and radii of the stellar 
components together with the individual apparent magnitudes 
determined from the light curve were used to estimate the dis- 
tance to DQ Vel (d = 3.112 kpc). 

Based on the position of the gainer on the HR diagram, which is 
well within the SPB instability strip, we interprete this periodic- 
ity of 5.14 d as a pulsation of an SPB type formed in the gainer. 
Similarities on the stellar parameters for DQ Vel and the DPV 
system V393 Sco were found during this study. However, we 
observed differences on the geometrical and physical properties 
of the two accretion discs, and we suggest that they are related 
to different evolutionary stages. We are currently working on a 
detailed analysis of the evolutionary stage of DQ Vel and aim to 
present our results in an forthcoming publication. 
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